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We present the results of a search for small-scale anisotropy in the distribution of arrival directions of cosmic 
rays above 10^^ eV measured in stereo by the High Resolution Fly's Eye (HiRes) experiment. Performing an 
autocorrelation scan in energy and angular separation, we find that the strongest correlation signal in the HiRes 
stereo data set recorded between December 1999 and January 2004 is consistent with the null hypothesis of 
isotropically distributed arrival directions. These results are compared to previous claims of significant small- 
scale clustering in the AGASA data set. 



1. Introduction 

Among the most striking astrophysical phe- 
nomena today is the existence of cosmic ray par- 
ticles with energies in excess of 10^° eV. While 
their presence has been confirmed by a number 
of experiments, it is not clear where and how 
these particles are accelerated to these energies 
and how they travel astronomical distances with- 
out substantial energy loss. 

At energies above several 10^^ eV, cosmic ray 
particles are believed to be of extragalactic ori- 
gin. If charged cosmic ray particles do not suf- 
fer considerable deflections in Galactic and ex- 
tragalactic magnetic fields, one can hope to iden- 
tify the sources and understand the underlying 
acceleration mechanism by a detailed study of 
their arrival directions. The strength and orienta- 
tion of these fields is poorly known and estimates 
vary |1I2| . but their impact should decrease at 
the largest energies; here, cosmic ray astronomy 
might be possible. 

The small cosmic ray data set has been sub- 
jected to extensive searches for an autocorrela- 
tion signal, i.e. clustering of arrival directions on 
small angular scales, and there have been a va- 
riety of attempts to correlate catalogs of known 
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astrophysical sources with cosmic ray arrival di- 
rections. 

So far, all efforts to identify the sources from 
the sparsely populated skymap have not pro- 
duced statistically convincing evidence for auto- 
correlation or correlations with any class of ob- 
jects. "Statistically convincing" should be em- 
phasized here, as there is actually no short- 
age of claims for both clustering and correlation 
with catalogs. Small-scale clustering of ultra- 
high energy cosmic rays above 4 x 10-*^^ eV has, 
for example, been repeatedly reported '3"4'5' 6|7] . 
with analyses mainly based on arrival directions 
recorded with the Akeno Giant Air Shower Ar- 
ray (AGASA) in Japan. If correct, these re- 
sults could indicate that cosmic rays originate in 
nearby, compact sources. 

However, there is considerable disagreement 
over the statistical significance of this cluster- 
ing signal. The problem arises from the way 
the chance probability of the signal is evaluated. 
Quite often, the data set used to formulate the 
correlation hypothesis is also used for evaluat- 
ing its significance. Problems with published 
claims of significant small-scale clustering have 
been pointed out by various authors |8l9ll0j . and 
it has become clear that ultimately, only statisti- 
cally independent data sets will allow a rigorous 
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Figure 1. Sidereal times of HiRes stereo events 
with energies E > 10 EeV. 



test of these claims. 

Such a data set has recently become available. 
The High Resolution Fly's Eye (HiRes) experi- 
ment in Utah has operated in stereo mode since 
December 1999. The operation of HiRes stereo 
and AGASA overlaps only a few months in time, 
but both detectors are located in the northern 
hemisphere and probe approximately the same 
part of the northern sky. The HiRes stereo data 
set therefore provides an opportunity to indepen- 
dently test previous claims that the arrival direc- 
tion of ultrahigh energy cosmic rays shows statis- 
tically significant small-scale clustering. 

2. HiRes Stereo 

HiRes is an air fluorescence experiment with 
two sites (HiRes 1 and 2) at the US Army Dugway 
Proving Ground in the Utah desert (112° W lon- 
gitude, 40° N latitude, vertical atmospheric depth 
860g/cm^). The two sites are separated by a dis- 
tance of 12.6 km. A general description of the 
HiRes detector is given in [TT] . 

Air fluorescence detectors observe air show- 
ers from the fluorescence light emitted by air 
molecules excited by particles in the shower. 
Tracks of air showers developing in the atmo- 
sphere are directly viewed by photomultiplier 
cameras that continuously watch the night sky. 

With showers detected out to distances of 
about 60 km, the instantaneous aperture of an 
air fluorescence telescope is quite large. However, 
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Figure 2. Angular resolution of HiRes stereo for 
energies 10, 40, and 100 EeV. The plot shows the 
fraction of simulated events reconstructed within S 
degrees of the true arrival direction, as a function 
of 5. 



the operating conditions of such an instrument - 
dark, moonless nights with a clear atmosphere - 
restrict its duty cycle to about 10%. Moreover, 
the sidereal time distribution (Fig.^l is modu- 
lated by overall seasonal variations as a result of 
the extended observation time during the winter 
months. 

In monocular reconstruction, the so-called 
shower detector plane, which is the plane that 
contains the shower track and the detector, is very 
well-defined, whereas the location of the shower 
within that plane has to be reconstructed from 
the trigger times of individual phototubes along 
the shower track. This leads to fairly large er- 
rors, of order several degrees, on the angle de- 
scribing the shower in the plane. The ambigu- 
ity is resolved if the shower is seen by two de- 
tectors. Stereo data has the distinct advantage 
that error bars on the arrival direction are fairly 
symmetric for most azimuth angles. A global y^- 



3 



10 EeV 



S 1 ^ 

o 0.8 - 

™ 0.6 ^ 

- 

^ 0.4 - 

O : 

0.2 - 



10 20 30 40 50 60 70 80 

zenith angle [°] 
1 EeV 



200 250 300 350 

azimuth angle [°] 



Figure 3. Angular resolution as a function of 
zenith angle (top) and azimuth angle (bottom) for 
simulated showers with energy E = IQEeV. The 
scatter plot shows individual events, the solid line 
represents the angular distance which includes 
68 % of all events in the zenith or azimuth an- 
gle bin. 



minimization using all available information, in- 
cluding the tube trigger times, gives an angular 
resolution which is typically about half a degree. 

The angular resolution of HiRes is determined 
using simulated showers. We use a full detec- 
tor simulation of proton showers generated with 
CORSIKA 6 [121 using QGSJET for the first in- 
teraction. The showers are thrown isotropically 
and undergo a full detector simulation and event 
reconstruction including all cuts that are applied 
to the real data As shown in Fig.|21 68% 

of all showers generated at 10^^ eV are recon- 
structed within less than 0.57° of the true shower 
direction. The angular resolution depends weakly 
on energy; the 68% error radius grows to 0.61° 
and 0.69° for showers generated at 4 • 10^^ eV and 
10^° eV, respectively, because at higher energies, 
showers are on average farther away. 



Fig. O shows the angular distance between true 
and reconstructed shower direction for a large 
number of showers with energy E = 10^® eV as a 
function of zenith and azimuth angle of the arrival 
direction. The angular resolution is essentially 
constant in zenith and azimuth angle, varying 
by less than 0.1° for zenith angles less than 70°. 
There is a small range of azimuthal angles where 
the angular reconstruction is notably worse: for 
azimuth angles where the shower and the two sites 
arc in the same plane, the stereo reconstruction 
looses its advantage and is essentially reduced to a 
monocular reconstruction. Showers having poor 
angular resolution as a result of this ambiguity 
fail the quality cuts. 

The major systematic errors on the arrival 
directions are due to uncertainties in the mir- 
ror pointing directions. Reconstruction of laser 
tracks and observation of UV bright stars in the 
field of view of the cameras indicate that this 
error is not larger than 0.2°. 

No explicit weather cut is applied for this data 
set. Simulations show that as long as weather 
conditions are known, their impact on the deter- 
mination of the shower geometry is minimal. The 
reconstruction uses an hourly atmospheric data 
base built from reconstructed laser shots. Note 
that stereo data provides us with a consistency 
check regarding the weather corrections applied 
to the data. If the distance of the shower to the 
two detectors is notably different, any inaccurate 
correction for light scattering and absorption will 
result in a mismatch of energy estimates for the 
two sites. 

After quality cuts the stereo data set 

used in this analysis comprises 271 events above 
10^^ eV recorded between December 1999 and 
January 2004. Fig.0] shows a skymap of their 
arrival directions in equatorial coordinates. 

3. Small-Scale Clustering 

As described in 13^, we search for small-scale 
clustering by performing an autocorrelation scan 
in energy and angular separation. For a given 
energy E^ we count the number of pairs sep- 
arated by less than 6. We then use sets of simu- 
lated showers with the same number of events as 
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Figure 4. Arrival directions in equatorial coordinates for 271 HiRes stereo events above 10 EeV recorded 
between December 1999 and January 2004- 



the HiRes stereo data set, but with an isotropic 
distribution of arrival directions, to estimate the 
chance probabihty P{E, 0) of finding Up or more 
pairs in a random data set just by chance. We 
identify the energy threshold Emm and the angu- 
lar distance dmin which gives the smallest chance 
probability, Pmm- This is the most promising 
potential clustering signal, and its statistical sig- 
nificance can be determined by performing the 
same scan over umc data sets with isotropic ar- 
rival directions, finding the minimum probability 
for each of these sets, and counting the number 
of data sets for which P < Pmin- Simulations 
show that despite the statistical penalty in- 
curred by scanning, this procedure results in a fi- 
nal chance probabihty of about 1 % (2 x 10^"^ %) 
for as few as 3 (4) clusters in 47 events, assuming 
that cosmic rays are not subject to strong deflec- 
tions in magnetic fields. 

For the HiRes stereo data set, we perform this 
scan over the 271 events with energies above 
10^^ eV. We scan from 0° to 5° in steps of 0.1° in 
angular separation. The energy threshold Emm is 



lowered one event at a time, starting at the high- 
est energy event and decrementing to 10^^ eV. 

For the HiRes stereo data set, the strongest po- 
tential signal is found at Emin = 1-69 x 10^^ eV 
and Omin = 2.2°. The chance probability for this 
potential signal is Pmm — 1-9 %, but the final 
chance probability after accounting for the scan 
is Pch = 52%. The signal is therefore not signifi- 
cant, and we conclude that the HiRes stereo data 
above 10^^ eV is consistent with the null hypoth- 
esis of isotropic arrival directions. As a conse- 
quence of our approach, we find that this conclu- 
sion is indeed very general. There is no evidence 
for significant clustering for any energy thresh- 
old above 10^^ eV on all angular scales of 5° or 
less. This indicates that clustering of arrival di- 
rections is weaker than previously suggested and 
might not be a general feature of ultrahigh energy 
cosmic rays after all. 

4. Comparison to the AGASA Data Set 

The HiRes stereo data set is still smaller than 
the published AGASA data set. In the energy 
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range of interest, above 4 x 10^^ eV, it has sig- 
nificantly less statistical power than the AGASA 
data set, but the difference in the number of 
events stems largely from a difference in observed 
flux rather than from a large difference in ex- 
posure. The discrepancy in flux is larger than 
the statistical errors of the experiments, and the 
possibility of a systematic energy shift of order 
30 % as studied in ^3] implies that by using the 
same energy threshold for both experiments, we 
are not comparing equivalent data sets. Our 
scan over energies is partly motivated by this en- 
ergy discrepancy. If the HiRes counterpart to 
the AGASA energy threshold is somewhat lower 
than 4 x 10^^ eV, then the scan at least ensures 
that we have tested for small-scale clustering in a 
HiRes data set equivalent in energy threshold to 
the published AGASA set. 

In order to evaluate the sensitivity of the cur- 
rent HiRes stereo data set to ultrahigh energy cos- 
mic ray point sources, we can set upper flux limits 
on sky locations which might harbor sources. The 
flux limits will be a function of the equatorial co- 
ordinates of any potential source, since the detec- 
tor exposure is not uniform. As an example, we 
determine upper flux limits for proton sources of 
constant intensity at the locations of the AGASA 
multiplets. 

To do this we track the "source" on its path 
over the (local) sky for every day of data tak- 
ing between December 1999 and January 2004, 
and simulate a constant rate of events with a dif- 
ferential energy spectrum c>c i?"^ ''^ from these 
positions using the correct zenith and azimuth 
angles. We use an hourly database of running 
times, weather conditions, and detector param- 
eters to generate events. After reconstruction, 
the quality cuts applied to the real data are also 
applied to the simulated events. We determine 
the number of events reconstructed within 2.5° of 
the "source" direction and compare this number 
to the generated flux to calculate a 90% upper 
flux limit using ^B] . TableQ] shows the result for 
the positions of the 5 AGASA multiplets above 
4 X 1019 eV. 

For AGASA, we calculate a flux from the loca- 
tions of the clusters by using the exposure pub- 
lished for the data set until May 2000. For each 
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HiRes stereo flux upper limits for the positions of 
the AGASA clusters. 
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AGASA cluster numbers according to ref. 0]. 



of the doublets (clusters 1, 3, 4, and 6 in Ta- 
bleUJ, the flux is about 0.010 km"^ yr~\ while 
for the triplet (cluster 2 in Table QJ, the flux 
is 0.013 km~^ yr~i. The HiRes stereo flux upper 
limits can at this point not rule out the flux es- 
timated from the AGASA clusters. Nevertheless, 
for some regions of the sky, the HiRes stereo expo- 
sure comes to within about 25 % of the exposure 
of the published AGASA data set. 

5. The AGASA Clustering Signal 

To understand the extent of the inconsistency 
between the HiRes result and AGASA's claim 
that ultrahigh energy cosmic ray arrival direc- 
tions exhibit significant small-scale clustering, we 
must better understand the AGASA clustering 
signal itself. An unbiased test of the clustering 
hypothesis is only possible if the data set used 
for formulating the hypothesis is not also used 
for evaluating the significance. Crucial parame- 
ters defining this signal, i.e. the energy threshold 
4 X 10^9 eV and the maximum angular distance 
between cluster events of 2.5°, were first described 
by the AGASA collaboration in Based on 
data recorded through October 1995, these cuts 
are justified a posteriori as cuts that lead to a 
significant clustering signal. Consequently, the 
strength of this signal should be evaluated using 
only data recorded since then. This split in "orig- 
inal" and "new" data was introduced in JOI , and 
the chance probability of the clustering signal in 
the new data set based on the cuts defined by the 
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original set was found to be 27%. If one allows 
for cross-correlations with the original set while 
replacing its clusters by single events at their av- 
eraged position, the statistical power of the new 
set can be increased, but the chance probability 
for clustering in the new set is still 8%, imply- 
ing that the clustering signal in the "new" set is 
consistent with the null hypothesis. 

One can take this independent test one step 
further. With 27 events, the new data set is 
slightly smaller than the original data set with 30 
events. If we extend the new data set by including 
the data taken since May 2000, summarized in a 
skymap on the AGASA web page ^17,, the origi- 
nal claim can now be tested with a larger data set 
comprising 42 events. The chance probability for 
the clustering observed in the enlarged new set is 
19%, so the larger data set fails to confirm the 
clustering hypothesis formulated using the origi- 
nal data set. 

These tests show that the persistent claims that 
the clustering signal has chance probabilities as 
low as 10~^ [1] are a direct consequence of the 
bias carried over from the original data set if the 
pre- and post-October 1995 data are not clearly 
separated. 

With evidence of small-scale clustering in the 
AGASA data set being weak to insignificant, our 
result is therefore consistent with AGASA's. It 
appears that at present, both data sets show no 
clustering signal. 

HiRes is currently the only cosmic ray detector 
probing the northern sky. The statistical power of 
the HiRes data set is increasing, and the impor- 
tant topic of small-scale anisotropics in the arrival 
directions of cosmic rays at the highest energies 
will be revisited once more data becomes avail- 
able. 
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